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^~>' The paper deals with Bianchi type V Universe, which has dynamical energy density. We consider Bianchi 

type V space-time, introducing three different skewness parameters along spatial directions to quantify the 
deviation of pressure from isotropy. To study the anisotropic nature of the dynamical dark energy, we 
assume that the skewness parameters are time dependent. It is found that the Universe achieves flatness in 
quintessence model. The physical behavior of the Universe has been discussed in detail. 
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• ^ I Recent observations have revolutionized our understanding of cosmology. Analysis of type la supernovae (SN la) 
(Perlmutter et al. 1997, 1998, 1999; Riess et al. 1998, 2004), cosmic microwave background (CMB) anisotropy 
<~i ' (Caldwell et al 2002; Huang et al. 2006), and large scale structure (Daniel et al. 2008) strongly indicate that 
C^. dark energy (DE) dominates the present Universe, causing cosmic acceleration. This acceleration is realized 
with negative pressure and positive energy density that violate the strong energy condition. This violation gives 
ps| . a reverse gravitational effect. Due to this effect, the Universe gets a jerk and the transition from the earlier 
' deceleration phase to the recent acceleration phase take place (Caldwell et al. 2006). The cause of this sudden 
On : transition and the source of accelerated expansion is still unknown. In physical cosmology and astronomy, the 
' simplest candidate for the DE is the cosmological constant (A), but it needs to be extremely fine-tuned to satisfy 
. the current value of the DE density, which is a serious problem. Alternatively, to explain the decay of the density, 
the different forms of dynamically changing DE with an effective equation of state (EoS) , lu = p^'^^' 
were proposed instead of the constant vacuum energy density. Other possible forms of DE include quintessence 
{uj > -1) (Steinhardt et al. 1999), phantom (w < -1) (Caldweh et al. 2002) etc. While the the possibility 
LU << —1 is ruled out by current cosmological data from SN la (Supernovae Legacy Survey, Gold sample of 
Hubble Space Telescope) (Riess et al. 2004; Astier et al. 2006), CMBR (WMAP, BOOMERANG) (Eisentein 
et al. 2005; MacTavish et al. 2006) and large scale structure (Sloan Digital Sky Survey) (Komatsu et al. 2009) 
data, the dynamically evolving DE crossing the phantom divide line (PDL) (w = — 1) is mildly favored. 

The anisotropy of the DE within the framework of Bianchi type space-times is found to be useful in gen- 
erating arbitrary ellipsoidality to the Universe, and to fine tune the observed CMBR anisotropics. Koivisto 
and Mota (2008a, 2008b) have investigated cosmological model with anisotropic EoS. They have proposed a 
different approach to resolve CMB anisotropy problem; even if the CMB formed isotropically at early time, it 
could be distorted by the direction dependent acceleration of the future Universe in such a way that it appears 
to us anomalous at the largest scales. They have investigated a cosmological model containing a DE compo- 
nent which has a non dynamical anisotropic EoS and interacts with the perfect fluid component. They have 
also suggested that cosmological models with anisotropic EoS can explain the quadrupole problem and can be 
tested by SN la data. Kumar and Singh (2010) have studied Bianchi type I cosmological models with constant 
deceleration parameter (DP) in the presence of anisotropic DE and perfect fluid. They have considered phe- 
nomenological parametrization of minimally interacting DE in terms of its EoS and time-dependent skewness 
parameters {S(t), 7(i), ?7(i)). Leon and Sarikadis (2010) have investigated that anisotropic geometries in modified 
gravitational frameworks present radically difference cosmological behaviors comparing to the simple isotropic 
scenarios. Akarsu et al. (2010) have investigated Bianchi-I anisotropic DE model with constant DP. Yadav and 
Yadav (2011a), Yadav et al. (2011b) have studied anisotropic DE models with variable EoS parameter. They 
have suggested that the dynamics of EoS parameter describe the present acceleration of Universe i.e. from earlier 
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deceleration phase to recent acceleration phase. Recently Pradhan et al (2011a, 2011b) have studied anisotropic 
DE models in different physical contexts. They have found that in the earlier stage EoS parameter was positive 
and it evolves with negative sign at present epoch. 

Bianchi type-V Universe is generalization of the open Universe in FRW cosmology and hence it's study is 
important in the study of DE models in Universe with non-zero curvature (1994). A number of authors such as 
Collins (1974), Maartens and Nel (1978), Wrainwright et al (1979), Canci et al (2001), Pradhan et al (2004), 
Singh et al. (2008), Yadav (2009) have studied Bianchi type-V model in different physical contexts. Recently 
Kumar and Yadav (2010) have studied isotropic DE model with variable EoS parameter in Bianchi type V space- 
time and found that the Universe is dominated by DE at present epoch and after dominance of DE, Universe 
achieves flatness. Following Eriksen et al (2004), it is found that some large-angle anomalies appear in CMB 
radiations which violate the statistical isotropy of the Universe. This motivates the researcher to consider the 
model of Universe with anisotropic DE. 

In this paper, we have studied some physically realistic and totally anisotropic Bianchi-V models with 
anisotropic DE and perfect fluid. To study the anisotropic nature of DE, we have assumed the time depen- 
dent skewness parameter, which modify EoS. The time dependent forms of the skewness parameter provide 
exact solutions of Einstein's field equation together with the special law of variation of Bubble's parameter. The 
paper is organized as follows. In Section 2, the models and field equations have been presented. The Section 3 
deals with the exact solutions of the field equations and physical behavior of the models. Finally, the results are 
discussed in section 4. 

2 Model and field equations 

The spatially homogeneous and anisotropic Bianca-V space-time is described by the line element 

ds^ = -dt^ + A^dx^ + e^°"={B^dy^ + C^dz^), (1) 

where A , B and C are the metric functions of cosmic time t and a is a constant. 

We define a — (ABC)^ as the average scale factor of the space-time (H]) so that the average Bubble's 
parameter reads as 
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where a — {ABC)^ is the average scale factor and an over dot denotes derivative with respect to the cosmic 
time t . 

The directional Bubble parameters along x , y and z coordinate axes, respectively, may be defined as 



Hx - ^, Hy - — , - —. (3) 



The Einstein's field equations (in gravitational units St:G = c = 1) read as 

R)-\9]R = -T^f"-T^f\ (4) 
where ' and Ty' ' are the energy momentum tensors of perfect fluid and DE, respectively. These are given 

T^f ' = diag (5) 



by 

T 

3 

and 



= diag [-1, ujx, ujy, uJz]p^'^'''' 

= diag [-1, w + (5, w + 7, w + 77]p('''=) (6) 

where p'^™-' and p*^™^ are, respectively the energy density and pressure of the perfect fluid component; p'^'*'^^ is 
the energy density of the DE component; S{t), j(t) and ri{t) are skewness parameters, which modify EoS (hence 
pressure) of the DE component and are functions of the cosmic time t; lu is the EoS parameter of DE; ujx, uiy 
and (jJz are the directional EoS parameters along x , y and z coordinate axes, respectively and we assume the 
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four velocity vector = (1, 0, 0, 0) satisfying u'^m = — 1. 

In a CO moving coordinate system (it' = Sq), the field equations for the anisotropic Bianchi type-I space- 
time (HI, in case of ([5]) and ([6]), read as 

+ (7) 
V")-(^ + 7)p(^^), (8) 
V™^ - (c^ + ry)p(^^\ (9) 

We assume that the perfect fluid and DE components interact minimally. Therefore, the energy momentum 
tensors of the two sources may be conserved separately. 

(™) J _ 
;j 

+ 3(p(™) +pi"^))H = 0, (12) 
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The energy conservation equation T^™' ■' — 0, of the perfect fluid leads to 



whereas the energy conservation equation t'^'^^^ ^ = Q, of the DE component yields 

■ide) ^ 3^(de)(^ ^ ^ ^(de) (^^^ ^ ^ ^^^^ ^ (^g) 

where we have used the equation of state p^''"^) = ujp'--'^'^^ . 

Equations ([7|)- ([TU| can be written in terms oi H, cr and q as 

pide) + 1 (3^ + J + ^ + ^) p{de) ^ jj2^2q - 1) - + ^ (14) 

Where q and a are deceleration parameter and shear scalar respectively. 

3 Solution of Field Equations 

Berman (1983), Berman and Gomide (1988), recently Kumar (2010) obtained some FRW cosmological models 
with constant DP and showed that the constant DP models stand adequately for our present view of different 
phases of the evolution of Universe. In this paper, we show how the constant deceleration parameter models 
with metric ([T]) behave in the presence of anisotropic DE. According to the law, the variation of the average 
Hubble parameter is given by Singh et al (2008) 

H = Da-"", (16) 

where D > and n > are constants. 



Following, Akarsu and Kilinc (2010), we split the conservation of energy momentum tensor of the DE into 

)ara 

p + 3p^'^'''> {lj + 1)H ^0 (17) 



two parts, One corresponds to deviations of EoS parameter and other is the deviation- free part of T^'^'^^ — 



and 

p^^'HdH,+^Hy+rjH,)^0, (18) 
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According to equations ([T7]) and ([T5)) the behaviour of p^^*^^ is controlled by the deviation-free part of EoS 
parameter of DE but deviations will affect p^'^'^^ indirectly, since, as can be seen later, they affect the value of 
EoS parameter. Of course, the choice of skewness parameters are quite arbitrary but, since we are looking for 
a physically viable models of Universe consistent with observations. We consider the skewness parameters i5, 
7 and 77 be function of cosmic time and we constrained (5, 7 and r\ by assuming a special dynamics which is 
consistent with eq. ([T5|). The dynamics of skewness parameters on x-axis, y-axis and z-axis are assumed to be 



-lit) = -aH^ 



Jde) ' 



(19) 
(20) 



r]{t) = -aHx- 



(de) 



(21) 



where a is an arbitrary constant, which parameterizes the anisotropy of the DE. In literature, many authors 
have considered totally anisotropic Bianchi-I, Bianchi-III and Bianchi-V space-times with only two skewness 
parameters of DE (Akarsu et al. 2010; Yadav et al. 2011a, 2011b; Pradhan et al. 2011a, 2Qllb). 

Finally, we assume that w = const, so that we can study different models related to the DE by choosing 
different values of w, viz. phantom (w < —1), cosmological constant (w = —1) and quintessence (cj > —1). 

In view of the assumptions ([T9 | - ([2T|) and uj ~ const., equation ([TT|) can be integrated to obtain 



(22) 



where po is a positive constant of integration. 

Integrating ([TT|) and absorbing the constant of integration in B or C, without loss of generality, we obtain 

= BC. (23) 

Subtracting ([7]) from ([8]), ([7]) from ([9]), ([8]) from (ITO|) and taking second integral of each, we get the following 
three relations respectively: 
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(24) 

(25) 
(26) 



where di , a;i , ^2 , 2^2 , and X3 are constants of integration. From equations ([Ml)- (1211) and (P^ . the metric 
functions can be explicitly written as 



A(t) — a exp 
B(t) — ma exp 



3uj 



I / a''-^dt + 



C{t) 



m ^aexp 



I / a^'^dt- 



where 



iuj 

(de 

3uj 

{X2 + X3) 



'+^^dt 



i-3(" + l)di 



with 

d2 = d^^, X2 = —xi. 

In the following subsections, we discuss the cosmologies for n ^ and ii — 0, respectively. 



(27) 
(28) 
(29) 

(30) 
(31) 



4 



3.1 DE Cosmology for n ^ 

In this case, integration of ([T6| leads to 



a(t) = {nDt)-, 

where the constant of integration has been omitted by assuming that a = at t = 0. 
Using (15^ into (P71) - (P^ . we get the foUowing expressions for scale factors: 



A{t) = {nDt)~ exp 
I 



2apo 



3(jjD(n — 3aj — 3) 



B{t) — m(nDt)" exp 
C{t) = m^^{nDt)" exp 



D{n - 3) 
-I 



{nDt)-r 
(nDt)- 



SuDln — 3aj — 3) 



{nDty 



-{nDty 



(32) 

(33) 

(34) 
(35) 



L»(n-3)' ' 3uD{n-Suj -3) 

where n 7^ 3. 

The physical parameters such as directional Hubble parameters {Hx, Hy, Hz), average Hubble parameter 
{H), anisotropy parameter (A), expansion scalar (6) and spatial volume (V) are, respectively, given by 



= {nty 



Hy = (nt)-i + l{nDt)^ + ^(ni^t)^ 
H = (nt)-\ 



-3(^ + 1) 



(36) 

(37) 

(38) 
(39) 
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V = {nDt)i exp(2Q;a;), 



Shear scalar of the model reads as 



= l\nDt) — + 



? 2 



3w^ 



(nDi)' 



The value of DP (g) is found to be 



aa 

q = -TT = n-l, 



(40) 

(41) 
(42) 

(43) 
(44) 



which is a constant. The sign of q indicates whether the model inflates or not. A positive sign of q, i.e., 
n > 1 corresponds to the standard decelerating model whereas the negative sign of q, i.e., < n < 1 indicates 
acceleration. The expansion of the Universe at a constant rate corresponds to n = 1, i.e., q = 0. Also, recent 
observations of SN la, reveal that the present Universe is accelerating and value of DP lies somewhere in the 
range — 1 < g < 0. It follows that in the derived model, one can choose the values of DP consistent with the 
observations. 

The skewness parameters of DE are as follows: 
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Figure 1: Plot of anisotropic parameter {A) versus time (t). 
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Figure 2: Plot of directional EoS parameters versus time (t). 



In view of ^ , the directional EoS parameters of DE are given by 
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The energy density and pressure of the DE components are obtained as 

p^^'^ = LopoinDt)^^ . 
From equations (fT4)) and p5|) , the pressure and energy density of the perfect fluid are obtained as 



(47) 
(48) 

(49) 
(50) 
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Figure 3: Single plot of energy conditions 
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The perfect fluid density parameter (ri*^™^) and DE density parameter {fl^'^'^'>) are given by 



(52) 
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Thus the overall density parameter (£7) is obtained as 



(53) 



(54) 



= 1 - 
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4q;Po 
3a;L)(n— 3a; — 3) 



(nDt)- 



.(55) 



It is observed that at t = 0, the spatial volume vanishes while all other parameters diverge. Thus the derived 
model starts expanding with big bang singularity at t = 0. This singularity is point type because the directional 
scale factors A{t), {B{t)) and C{t) vanish at initial moment. The solutions for the scale factors have a combina- 
tion of a power-law term and exponential term in the product form. The DE term appears in exponential form 
and thus affects their evolution significantly. For a < and n > 3u; -I- 3, the DE contributes to the expansion of 
A{t) while opposing to the expansion of B{t) and C{t). Likewise, for a > and n > 3a; -I- 3, the anisotropic DE 
opposes the expansion of A{t) while contributing to the expansion of B{t) and C{t). 

The difference between the directional EoS parameters and hence the pressures of the DE, along x-axis and 
y-axis (or z-axis) is 3a{nt)~^, which decreases as t increases. Therefore, the anisotropy of the DE decreases as 
t increases and finally drops to zero at late time. The variation of mean anisotropic parameter (A) versus has 
been graphed in Fig.l by choosing D = 2,uj = —l.l,n — 0.5 and other constant as unity. Since the current 
observations strongly recommend that the present Universe is accelerating (i. e. q < ). We consider n = 0.5 
i. e. g = —0.5 in the remaining discussion of the model. 

Fig. 2 depicts the variation of directional EoS parameter (i.e. uix^ ujy, luz) versus cosmic time. We observe 
that directional EoS parameter along x-axis {i.e. ujx) is decreasing function of time while directional EoS param- 
eters along y-axis (or z-axis) are increasing function of time. At the later stage of evolution, all the directional 



7 



-20 - 
-26 - 
-30 - 
-35 
-40 
-45 



!)l"'> 



Figure 4: Plot of density parameters versus time (t). 



EoS parameters approaches to —1 as expected. The same is predicted by current observations. 

FoUowing, Caldweh (2002), Srivastava (2005) and recently Yadav (2011) have investigated phantom model 
with uj < —1. They have suggested that at late time, phantom energy has appeared as a potential DE candidate 
which violates the weak as well as strong energy condition. The left hand side of energy conditions have graphed 
in Fig. 3. 

From Fig. 3, for w = —0.9 (i. e. quintessence model), we observe that 

(i) p^'^"^ > 

(ii) -f- > 

(iii) p^'^"') - > 

(iv) p^''^) + 3p(''^) < 

Thus the derived quintessence model violates the strong energy conditions, as expected. 
Further, for uj — —l.l (i. e. phantom model), it is observed that 

(i) p^'^"^ > 

(ii) - p^"*^) > 

(iii) +p(d'^) < 

(iv) p^**^) + Sp^''^) < 

Thus the derived phantom model violates the weak as well as strong energy conditions. The same is predicted 
by current astronomical observations. 

From equation ([55|) . it is observed that for < n < 1 and uj > —1, the overall density parameter (ft) 
approaches to 1 for sufficiently large times. Thus, the derived model predicts a flat Universe in quintessence 
model for sufficiently large times. Fig. 4, depicts the variation of density parameters versus cosmic time during 
the evolution of Universe. 
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3.2 DE Cosmology for n = 

In this case, integration of ([T6| yields 



a{t) — cie 



Dt 



where Ci is a positive constant of integration. 
The metric functions, therefore, read as 
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The other cosmological parameters of the model have the following expressions: 



3(w+l) 



Hy = D 



apo -3D(uj+l)t 



3(w+l) 



H -PI "P" -3D(u; + l)t 
^^^^+0 3(.^+l)^ 

SuJCi 
H = D, 

2a^Po' 



A = 



o 7-i2 2 6(ti;+l) 



-6D{ui+l)t 



e = 3D, 

^2 ^ Q^^PO^ ^-6D(uj+l)t 
„ 9 6(w + l) 

3w^Cj 

The skewness parameters and the directional EoS parameters of DE, respectively, are given by 
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Figure 5: Plot of directional EoS parameters versus time (t). 
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Figure 6: Plot of Anisotropic parameter [A) versus time (t). for n=0 
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Figure 7: Single plot of energy conditions in quintessence model for n = 
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Figure 8: Single plot of energy conditions in phantom model for n = 0. 
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The energy density and pressure of the DE components are given by 



(74) 



(75) 
(76) 



From equations (fT4|) and (fT5|) . the pressure and energy density of the perfect fluid components are obtained as 
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The perfect fluid density parameter (£7*^™^) and DE density parameter {fl^'^'^'>) are given by 
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Thus the overall density parameter (Jl) is obtained as 



n = 1 
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a Po 



(77) 
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We observe that the model has no initial singularity. The directional scale factors and all other physical 
quantities are constants at i = 0. The directional scale factors and volume of the universe increase exponentially 
with the cosmic time whereas the mean Hubble parameter and expansion scalar are constants throughout the 
evolution. Therefore, uniform exponential expansion takes place. Further, we see that the DE term appears in 
exponential form in the scale factors and thus effects their evolution significantly. Thus, the spatial geometry of 
the universe is affected by the anisotropic DE. The difference between the directional EoS parameters of the DE 
and hence the pressures of the DE, along x-axis and y-axis (or z-axis) is SaD, which is constant throughout the 
evolution of the universe. Therefore, the anisotropy of the DE does not vanish during the evolution of universe. 

As t ^ oo, the scale factors and volume of the universe become infinitely large whereas the skewness and 
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directional EoS parameters, directional Hubble parameters become constants. The pressure and energy density 
of the DE drops to zero in quintessence model (i. e. w > —1). Fig. 6 shows the plot of anisotropy param- 
eter (A) versus time, in both quintessence model and phantom model. We see that in quintessence model, 
the anisotropy parameter decreases as time increases and finally drops to zero at late time. But in phantom 
model the anisotropy parameter does not vanish during the evolution of Universe. The left hand side of energy 
conditions have graphed in Fig. 7 and Fig. 8 for quintessence model and phantom model respectively. It is 
observed that the quintessence model violates strong energy condition whereas phantom model violate weak 
energy condition as well as strong energy condition, as expected. 

For n = 0, we get q = —1; incidentally this value of deceleration parameter leads to ^ = 0, which 
implies the greatest value of Hubble's parameter and the fastest rate of expansion of the universe. Therefore, 
the derived model can be utilized to describe the dynamics of the late time evolution of the actual Universe. So, 
in what follows, we emphasize upon the late time behaviour of the derived model. Fig. 5 depicts the variation 
of directional EoS parameter versus time for n = 0. We observe that uJx,i^y or lOz are evolving with negative 
sign, as expected. 

From equation (|5T|) . it is observed that for sufficiently large time, the overall density parameter ($7) ap- 
proaches to 1 in quintessence model. Thus the derived model predicts a flat Universe at late time. 



4 Concluding Remarks 

In this paper, some spatially homogeneous and anisotropic DE models in Bianchi type V space-time have been 
studied. The main features of the work are as follows 

• The models are based on the exact solution of Einstein's field equations for anisotropic Bianchi type V 
space-time filled with DE. 

• The singular model {n ^ 0) seems to describe the dynamics of the Universe from big bang to the present 
epoch while non singular model (n — 0) seems reasonable to project dynamics of the future Universe. 

• The directional EoS parameters (i. e. uj^, ujy or uJz) evolve with in the range predicted by observations. 

• In the present models, we do not rule out the anisotropic nature of DE. The anisotropic DE contributes to 
the expansion of one (or two) of the scale factors while it opposes the expansion of the other two (or one) 
scale factors leading the geometry of Universe. The anisotropy of DE vanishes at late time for singular 
model (see Fig. 1) whereas for non singular model, anisotropy occurs at early stage (i. e. quintessence 
model) or at later time of Universe (i. e. phantom model) (see Fig. 6). 

• The derived quintessence models violate the strong energy condition whereas the phantom models violate 
the weak energy condition as well as strong energy condition (see Fig. 3, Fig. 7, Fig. 8). 

• The flatness of Universe can be achieved in quintessence model {i.e. uj > ~1) for sufficiently large time 
because the overall density parameter (f2) approaches to 1. Thus in our analysis, the quintessence model 
is turning out as a suitable model for describing the late time evolution of Universe. 
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